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Abstract 

We discuss a particular secluded WIMP dark matter model consisting of neutral 
fermions as the dark matter candidate and a Proca-Wentzel (PW) field as a media¬ 
tor. In the model that we consider here, dark matter WIMPs interact with standard 
model (SM) particles only through the PW field of ~ MeV - multi-GeV mass par¬ 
ticles. The interactions occur via a U{1)' mediator, V^, which couples to the SM 
by kinetic mixing with U{\) hypercharge bosons, One important difference 
between our model and other such models in the literature is the absence of an ex¬ 
tra singlet scalar, so that the parameter with dimension of mass My is not related 
to a spontaneous symmetry breaking. This QED based model is also renormal- 
izable. The mass scale of the mediator and the absence of the singlet scalar can 
lead to interesting astrophysical signatures. The dominant annihilation channels 
are different from those usually considered in previous work. We show that the 
GeV-energy y-ray excess in the galactic center region, as derived from Fermi-hAI 
Gamma-ray Space Telescope data, can be attributed to such secluded dark matter 
WIMPs, given parameters of the model that are consistent with both the cosmo¬ 
logical dark matter density and the upper limits on WIMP spin-independent elas¬ 
tic scattering. Secluded WIMP models are also consistent with suggested upper 
limits on a DM contribution to the cosmic-ray antiproton flux. 
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1. Introduction 

The clear astronomical [[ll] and cosmological evidence for large amounts of 
dark matter (DM) in the universe [Q] has led to the construction of various theo¬ 
retical models that go beyond the standard model (SM) weak-scale theories and 
which attempt to account for the DM abundance in the universe . The observa¬ 
tional evidence for DM has motivated various experimental searches to find dark 
matter [0]. 

Data from colliders are used to search for evidence of dark matter particles. 
Experiments with detectors like DAMA, CoGeNT, CDMS, XENON and FUX 
are used to search for evidence of the recoil en ergy of nuclei that would be pro¬ 
duced by scattering with dark matter particles O, |^. High-energy colliders like 
FHC (Farge Hadron Collider), have obtained significant upper limits on the anni¬ 
hilation of WIMPs to quarks [13] • They also offer very interesting possibilities to 
investigate interactions involving DM mediators. 

Space-borne detectors have been used to search for evidence of the products of 
of dark matter annihilation, particularly y-rays and cosmic-ray positrons. These 
searches have conservatively produced constraints on dark matter annihilation, 
both from cosmic y-ray studies Jsl] and cosmic-ray positron studies [0] ■ However, 
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analyses of the Fermi-LP^ data have indieated the existenee of an ’’excess” flux of 
y-rays above that expected from cosmic rays interacting with interstellar gas. This 
flux appears to be extended around the region of the galactic center. It appears to 
peak in the 2-3 GeV energy range. This excess has been interpreted to be a 
possible indication of the annihilation of weakly interacting dark matter (WIMP) 
particles having a mass in the 20 - 45 GeV range, annihilating primarily into 
quarks or, less likely, in the 7-12 GeV mass range annihilating primarily into 
charged t leptons fiaini] . Models involving other annihilation channels have also 
been considered [|l2ll . We note, however, the determination of various possible 
components of y-ray emission from the galactic center is complicated and the 
y-ray data are not precise enough to point to a unique origin. Other possible 
interpretations of this excess that include other contributions to the y-ray flux in 
the region of the galactic center have been suggested Q. 

Neutralino supersymmetric WIMPs, viz., the lightest supersymmetric dark 
matter particles, have been a popular choice to be the DM WIMPs because they 
are stable and neutral and their cross section naturally leads to the correct cosmo¬ 
logical DM density. However, as of now, the LHC has not found any evidence for 
such particles. Therefore, other candidate WIMP models have been explored and 
should be further explored. 

In this work we consider a model of WIMP dark matter in which the dark sec¬ 
tor is just quantum electrodynamics (QED) extended with a new massive photon 
field, usually dubbed a Proca-Wentzel (PW) field. It is well known that this is a 
renormalizable theory because it couples to a conserved vector current. Hence, 
the dark sector is made up of Dirac fermions, rj, that only interact via a PW field, 
here denoted by V'^, that serves to mediate between DM fermions and standard 
model particles. 

Diagonalization of the kinetic terms in the Lagrangian give both the standard 
Z boson and an extra neutral gauge boson that we denote as Z'. We will show that 
this DM model can produce the observed cosmological DM abundance and that 
its annihilation into standard model (SM) particles can also lead to astronomically 
observable fluxes of y-rays in the galactic center and in dwarf galaxies. 

In our model secluded DM interactions occur only through Z' mediators which 
subsequently decay to SM particles. Thus, they have an very small elastic scat¬ 
tering cross section with nuclei. This distinguishes such secluded WIMP models 
from other WIMP models. The decay of Z'’s into quark-antiquark channels pro¬ 
duces pions, among which are 7r°’s that decay to produce y-rays. The .7r°-decay 


) proi 

y-ray spectrum has a characteristic peak at m„/2 [Il4ll and is bounded by the rest- 
mass of the WIMP, m,i (e.g,, llisll ). The Z' decay, particularly into charged leptons 
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and light quarks, also yields y-rays through internal bremsstrahlung. In this pro- 
eess, a y-ray speetrum is produeed that peaks near [Il6ll . Eleetrons resulting 
from this proeess ean produee y-rays via Compton seattering in the interstellar 
medium. 

Other Z' models have reeently been diseussed in the eontext of astrophysieal 
y-ray produetion (see, e.g.. Ref. lll7ll '). However, our model, in whieh the Z' is the 
mediator of the DM interaetions through the annihilation ehannel rj + 1} ^ Z' + Z', 


as shown in Figure 1, and with niz' m,,, was not eonsidered in Ref. [Il7ll . (See 
also the diseussion in SeetionH]). 

The outline of the paper is as follows: In seetion |2] we present the model. 
In seetion |4] we diseuss the differenees between our model and previously ex¬ 
plored WIMP DM models. In seetionwe ealeulate the relie density of our Dirae 
fermion DM eandidate. In seetion |7] we diseuss the astrophysieal produetion of 
y-rays from from annihilation of our seeluded WIMPs as a possible explanation 
of the so-ealled y-ray exeess from the galaetie eenter and potential y-ray signals 
from the Milky Way satellite dwarf galaxies. In seetion [8] we summarize our re¬ 


sults and eonelusions. In appendix Appendix A we give the full eouplings for the 
interaetions deseribed by the model. 


2. The Basic Model 

The possible existenee of an extra [/(!)' symmetry of nature beyond the SM 
has been eonsidered for a long time. The addition of this new symmetry fae- 
tor to the eleetroweak 5 C/(2) (g> f/(l) of the SM oeeurs via the so-ealled kinetie 
mixing portal, by mixing with the hypereharge gauge boson A real massive 
veetor field eoupled to a fermionie veetor eurrent is a well behaved theory. The 
mass term breaks what would-be a gauge symmetry (whieh is valid in the kinetic 
terrn). However, this only implies the constraint = 0. For earlier references 
see [18]. 

The dark matter fermion WIMPs of the model, rj, interact only with the PW 
field that is the connection between DM and ordinary matter, designated by V'^. 
The PW field mixes through the kinetic term with the U{1)y vector field of the 
SM, at this stage designated by B'^. DM interactions occur only through Z' medi¬ 
ators which subsequently decay to SM particles. As such, secluded DM WIMPs 
have a very small elastic scattering cross section with nuclei. This distinguishes 
secluded WIMP models from other WIMP models, thus allowing for a potential 
experimental test. 
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After the diagonalization of the mixing in the kinetie terms, 


V' 

M 




s\b 






( 1 ) 


where is a linear eombination of Z and the extra neutral gauge boson Z'. Cos¬ 
mological DM abundance constraints on the DM annihilation cross section will 
then require a small mixing angle between V and Z' so that V - Z' [[l9ll . Depend¬ 
ing on the mass of Z', interesting signatures for these types of mediators could 
come from the Drell-Yan channel pp —> Z' —> ZZ and from non-conventional de¬ 
cays of SM Higgs boson such as /? —> ZZ', potentially observable with the Large 
Hadron Collider (LHC). In this work we consider a scenario where the Z' is light, 
such that Mz' ^ The scalar sector is the same as that of the SM, viz., a doublet 
with Y = -1-1. 

Empirical bounds on Z' couplings exist for these kinds of models. Such 
bounds depend on the mass scale of Z'. High energy colliders are sensitive to 


Mz 


> 


10 GeV and the constraints for lighter Z'’s are given mainly by preci 


sion QED observables, B meson decay and some fixed target experiments 112ih . 
Besides these constraints, there is also a constraint on the lifetime of a Z'. Its life¬ 
time should be less than one second in order to guarantee that the Z' decays before 
the onset of big-bang nucleosynthesis ll22ll . As a consequence, the secluded DM 
fermion will annihilate preferably into a Z' pair 
The full Eagrangian of the model is given by 


~ J-^SM J~Dark -i^Dark+mt^ ( 2 ) 

where £,sm is the SM Eagrangian, and Jlnark is the dark Eagrangian given by 

■Loark = - nir,)!!. (3) 

The mixing between dark and SM matter, occurs through the third term in Eq. Q, 
which is 


-C'Dark^int = g.^rriV'^ + ^MlV'^V'^ (4) 

- — V' A- 

4 2 4 ’ 

where with = V'^, B'^, and B'^ is the abelian gauge boson of 

the U{1)y factor in the SM Eagrangian. In Eqs. ([3]) and dH), the masses m,j,Mv 
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and the couplings gvB and gr, are free parameters. The mixing in the kinetic term 
in Eq. (01) can be diagonalized by a GL(2, R) transformation [l23n 


V 

B 




SvB 
SVB 1 


y' 

B' 


(5) 


Using ([5]) in (01) we obtain 


£! 


Dark+int — 


Ti'fr]Vf, + 


Ml 

^ ■y„y^ 


2(l-gt'a) 


2 


( 6 ) 


yj^-slB 

After symmetry breaking we get the mass matrix in the (5^, and V^) basis 




2 2 

g Vh 

4c^ 




0 

0 

0 


0 0 
1 -^Sw 


U, 


(7) 


where sw is the usual weak mixing angle, ^ = gvB/ V> - gig, r = Ml!Ml, and 
where the matrix U is given by 


U 


C\f/ s^/|/ 0 ^ 

— Syv Cyv 0 

0 0 1 


We can write the 3x3 mass eigenstates matrix as 


( 8 ) 


with 


' B ' 


Cw 

^W^a 


\ 

' A ' 

lys 


Sw 

CwCa 

~Cw 


Z 

1 J 


[ 0 

Sa 


J 

7' 

\ ^ ) 


ha - -- 


2sw^ 


Vs 


( 10 ) 


where ha denotes tan 26^, U denotes tan 6^, denotes cos 6a and Sa denotes the 
sin 9a. For small values of da, we can expand ta - hal2.-t\^l'& and use the relation 
Cff = (1 + to perform the calculations. 
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The masses of the Z and Z' are found by diagonalizing the matrix (|7]). They 
are given by: 

Ml ^ ^ I - 

Mz,z' = ^[(1 + + ^) ± V(i - (11) 

where Mzg = gVhl(2cw), assuming that ^ 1 and r < (1 - (Again, 

we assume that the lighter neutral veetor boson is Z'). Notiee that for the elee- 
troweak preeision observables, it follows that as we inerease Mz' we have to 
deerease the eoupling gys in order to respeet the eonstraint whieh implies that 
< 10 - 2 . 


3. Classification of Heavy Photon Models 

Eleetroweak models with an extra t/(l) symmetry are among most well moti¬ 
vated extensions of the Standard Model. In a general eontext we use the notation 
f/(l)i( 8 )f/(l) 2 . All these models allow a kinetic mixing between the field strength 
tensors of both [/(I) gauge bosons, F\^yF^. 

We can separate this sort of models in several groups. For instance, 

A1 Both [/(I) groups are visables. It means that SM fermions carry both quan¬ 
tum numbers. This possibility usually arises in the context of grand unified 
theories (GUTs). Examples are the models of Babu et al. and that of Galison 
and Manohar ifl^ . Models in which U{\)\ = 1/(1)^ and U{1)2 = U{1)b-l 
are of this type. 

A2 There are fermions which carry only one of the t/(l) charges, others carry 
both charges and some with no charges. As an example of these sort of 
models we have del Aguila et al. in [jisl] . 

A3 One of the f/(l) factor is visible, the SM particles carry one of the f/(l) 
charge, for example f/(l)i, an the other, [ 7 ( 1)2 is dark. Our model is of this 
type. 


Another way to classify these sort of models is by considering the mass scale 
at which the kinetic mixing occurs. Although it is not usually explicitly say, this 
is an important point. 


B1 Models in which [/(l)i = U{1)y. The kinetic mixing occurs before the SSB. 
Our model and the one of Ref. 112311 are of this type. 


8 






B2 Models in which U{\)2 = U{\)q, i.e., the kinetic mixing occurs after the 
SSB and the mixing is with the massless photon. An example of this type 
of model is that of Holdom in lllSn . 


B3 Neither t/(l)i nor [7(1)2 are related to the symmetry of the SM directly, but 
[/(l)i <S) [7(1)2 —^ U{\)y after SSB. The Galison and Manohar model |i8] is 
of this type. 

There are also models in which 

Cl The mixing among the two [7(1) factor occur in the ki netic and also in 
the mass terms. In this case the photon has a component on Z'. See, for 
instance, the model of Babu et ah, in Ref. lllSn . 


C2 The mixing among the two [7(1) factors occurs only in the kinetic term. 
Examples of these sort of model are our model and that of Ref. ll23ll . In 
these models the photon has no component on Z'. 

We can also classify the models according to 

D1 Extra scalars, for instance singlets, are added to break the additional [7(1) 
symmetry. 

D2 The only scalar in the model is that of the SM. Our model is of this type. 

Models may also be classified according to the fermion, scalar or vector nature of 
both the DM and the mediator. 

El Our model has a Dirac fermion as DM and a real massive vector as mediator. 


The model of Ref. [|24ll postulates a real scalar or a vector field to be the dark 
force and a complex scalar or a fermion to be the dark matter. 

El The Stiickelberg scalar model was introduced in order to maintain the gauge 
invariance in a QED-like theory with a massive real vector field. In this 
sense it provides an alternative to the Higgs mechanism. However, it needs 
an axionic scalar S in the notation of Ref. KSIB . 
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4. Specific differences between other modeis and our modei 


Our model is motivated by the following eonsideration. It is well known that 
massive quantum eleetrodynamie is a renormalizable and free of anomalies theory. 
It means that spontaneous symmetry breaking in order to give mass to the veetor 
field ean be implemented, but it is not mandatory. Our model has a Dirae fermion 
as DM and a real massive veetor as mediator. 

The important point is that the photon (massless or not) eouples to a eonserved 
eurrent. This oeeurs beeause the bad high energy behavior of the veetor propaga¬ 
tor ~ kf^kylniy vanishes when eontraeted with the eonserved eurrent, say J^. The 
U(iy symmetry in the kinetics term is broken by the mass term (my/2)V^y^, but 
its only consequence is the constraint df^V^ = 0, which comes from the equation 
of motion. It is still possible to restore the gauge invariance if a massive scalar 
field called a Stiickelberg field is invoked. If S {x) denotes the gauge Stiickelberg 
scalar, satisfying the equation (□ - m\)S {x) = 0. 

We note that the Proca-Wentzel Lagrangian that we use is, in fact, equivalent 
to the Stiickelberg Lagrangian when S {x) = 0 (For the main references, see [H^). 
Thus, our model coincides with that of Ref. [1251] (FI above) when no Stiickelberg 
field is present. An important difference between our model and a Stiickelberg- 
like model is that the later predicts the existence of millicharged particles [|27|] . 
That model introduces a new unit of electric charge, allowing the electric charge 
to have non-integer values. Hence, although both models have similar features 
there are many differences between the models. 

In our model we have a Dirac fermion and a vector mediator and the mediator 
is lighter than the DM candidate. This is also the case in Ref. lE^ . However, 
in the later model, the leading interaction between the Standard Model and the 
dark sector is the kinetic mixing between the photon and the U'{1) gauge boson: 
X = (l/2)eF))yF^’'. The leading kinetic mixing with the photon is also used in 
Ref. m- 


5 . Details of the model 

In the present model, the kinetic mixing is with the gauge boson of the U{1)y 
factor before the spontaneous symmetry breaking induced by the Higgs boson in 
the SM that we have denoted B'. It is in this state that the mixing in Eq. (5) 
occurs. Hence, the leading interaction between the SM particles and the dark 
sector is through the and B component on Z'. See equation (9). It is through 
the mixing with the Z boson that the particles in the SM have the couplings with 
Z' given in the Appendix. 
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A model in which the photon has no dark component is that in Ref. [12311 . 
however, they introduce an scalar singlet O//, to break the f/(l)' spontaneously. 
We denote the vacuum expectation value of this singlet vs, as opposed to the VEV 
of the SM Higgs, which we denote by Vh . The scalar singlet make the difference 
of the model in Ref. [|23n and ours. In former case the interactions the of Higgs 
bosons are 

nif 

hff : -ich—, 

Vh 
^2 


hWW : 2m 


m 


h- 


IV 


Vh 

m| niy 

HTjTj . 'XiCji 

Vh Vs 

ml My 

hZ Z : 2iCh—{Sa + ^swc„) - iSh —c„, 

Vh Vs 

2 2 

rriy 

hZZ . 2icji [(Cq. isii s^Cq. 


Vh 


Vs 


where Vs is the VEV of the scalar singlet (designated by ^ in Ref. I123h l that is 
needed to break the t/(l)' symmetry and cu = cos 0/„ s/, = sin0/, are the cosine and 
sine of the mixing angle: 


/ (psM ^ \ / /* \ 

\ 4>h ) \ -Sh Ch ) \ H ) 


( 12 ) 


with h and H being the mass eigenstates. The scalar singlet carries U{\y charge 
denoted e' in [Il9n . 

Moreover, the scalar potential has several dimensionless parameters in the 
quartic interaction terms 


y(^)( 0 , 5 ) = d(Oj^OsM)' +p(S^ 5)2 + /cS^Oj^OsM 


which obeys the relation [|23h . 

tan(20/,) 


/rVhVs 


- dVh^ 


(13) 


(14) 


In a model with a scalar singlet, the SM-like Higgs scalar has weaker inter¬ 
actions with fermion and vector bosons. Such interactions are suppressed by a 
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1 


1 


rwwz' 

'n 

WWZ' 


Figure 1: Feynman diagram for rj + rj ^ Z' + Z'. 


factor of Ch relative to SM interaetions. The vertiees hZZ, hZZ' and hZ'Z' have an 
extra term proportional to s^. Only when c/, = 1, do the interactions with the SM 
partieles for both models agree. However, in the models of Refs, [fl^ and [231. 
the dark veetor still interacts with the singlet h. In general, the phenomenology 
of both models is a bit different in proeesses like h —> Z'Z, Z'Z',ZZ [21]. This 
may imply that the region allowed for the four parameters mrj,g,j,gvB and my is 
different in both models, although an overlap may also exist, as follows from the 
DM ealeulation above. For example, in the models of Refs. ifl^ and i^, the 
mass of the veetor field my = g^^Vj and the mass of h also depends on V 5 . Henee, 
eonstraints on my are also eonstraint on m/, which could be in disagreement with 
aeeelerator data. This relation does not exist in our model in whieh there is no 
sealar singlet. In our model my is a free parameter. These interaetions are equal 
to those in the appendix only when c/, = 1 in equation (IT^ . 

Notiee that we are eonsidering the Z' to be lighter than the Z, then in our 
Eq. (11) the signal - eorresponds to the lighter boson. In faet, it may have zero 
mass if My = 0. In this case we have a dark QED indeed. We also note that the 
values that we take for gvB (called e in experimental papers) are mueh lower than 
the present experimental upper limits Hi but may be within the sensitivity range 
of future experiments [l32n . 


6. Cosmological DM Density 

As noted before, in this model, in addition to SM partieles, there is a dark 
seetor eomposed of a Dirac fermion, 77 , and a Proea-Wentzel veetor field, - Z', 
in the limit of small mixing angle. We eonsider rj to be the only eomponent of 
DM. In the regime where m,, » Mz', the primary annihilation channel will be 
T] + 7 ] ^ Z' + Z' [Il9n . as shown in Figured] The differential annihilation eross 


12 










section in this case is given by 

do-gnn _ I S 

dO.' 64;r5(l - s - 

[Ai + cos 9[-AI + AI cos 6(Am^ + s + A 2 cos 6)]] 
5(^2 cos 6 - s)'^{s + A 2 cos 0) 


( 15 ) 


where dQ.' = dcp sin 6d6, Ai = 5(32m^ - Sm^s - s^) and A 2 = 


^s{s - 4m2), 


where the summation over final spins and average over initial spins are taken into 
account. 

After integrating equation (fTSl) to obtain the total annihilation cross section, 
we used the approximation for the square of the center-of-momentum energy, s - 
4m^ + m^v^, which is valid for non-relativistic particles. The annihilation cross- 
section can then be expressed in the form (cTann I v |) = a -l- bv^, where the a and 
b are given in equation (fT^ . Note that relation between the velocity in the center- 
of-momentum (Vcms) and the relative velocity (v) is given by Vans = v/2 [H^. We 
then find 


ann^'} 


4 4 
Cagn 


3 

+ -■ 


4 4 
Cagn 


(16) 


32;r(l - gy^W^ 8 32;r(l - gl^W^ 

Using this expression, we calculated relic density as a numerical solution to the 
Boltzmann equation, discussed in [134 13511 




1.04X lO^GeV-' X 


/ 


M, 


pi 


where Xf is given by 


Af = In 


c(c + 2) 


^/^JXfja + 3blXf 


45 gm^Mpiia + bbjXf) 


(17) 



(18) 


and where g = 2 for fermionic DM, the Planck max, Mpi = 1.22 x 10^^ GeV, 
c is a parameter of order unity considered here as 5/4 and g^Q^f) is the number 
of relativistic degrees of freedom at freeze-out. Given one neutral Dirac massive 
particle and one additional neutral gauge boson and all the SM content, g»(A/) « 
113.25, Xf = nirilTf and Tf is the temperature at freeze-out. For relics with mass 
in the range of electroweak scale. A/ is in the range 20 - 30. 

Another parameterization for {(TannV), given in Refs. S and ll^ . is often 
used to obtain the annihilation cross section for fermionic dark matter. In the case 
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Figure 2; {(Tamv) in crn^Is as a function of m,j in the limit where Mz’ Im^ —> 0, ^ ~ gvB = 8x10 

g, = 0.088. 

where the mass of the mediator is negligible, sinee gvB ^ 1, our equation (fT^ 
reduees to equation (6) in Ref. viz., {(TannV) = {I/2)7r{a'^/ It follows 
from the identifieation of these two equations that the value of the parameter that 
determines the eorreet DM abundanee, viz. a', is aetually independent of m,,. Con¬ 
sidering that the annihilation in the early universe is predominantly s-wave, we ean 
negleet the b term in equation (fT^ in the s-wave approximation and identify 



(19) 


a 


Figures [2]-[5] give the results for < cra„„v > and Q.h^ as a funetion of in the 
limit where Mz' /m^j 0, taking the values for the model parameters as indieated 
in the figure eaptions. Figure [6] shows the relation between and whieh yields 
the eorreet DM abundanee to within 3cr. The spread in the eurves is determined 
by taking Q.jjm = 0.1196 and a Scr spread with cr = 0.0031. 

7. Gamma-rays from t] Annihilation 

We have pointed out that our partieular seeluded WIMP dark matter model 
has eertain attraetive theoretieal features. Having further shown that our proposed 
DM eandidate ean yield the proper eosmologieal relie density, we next eonsider 
the produetion of astrophysieal y-rays from rj annihilation. 


'The factor 1/2 in equation this equation comes from averaging over the t/(l)' charges. 
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Figure 3; Q.h^ as a function of in the limit where Mz’ Iniri ^ 0, ^ ~ gvB = 8x10 gjj - 0.088. 
The horizontal lines denote the range of Q.h^ measured by Planck Collaboration 11361] . 



Figure 4; {cTannv} in cm^ j s as a function of ^ ~ gvB = 10'’, ^^, = 1.12x10 Mz- =3.5 GeV. 



Figure 5 : Q.h^ as a function of ^ ~ gvB = 10“®, grj = 1.12 x 10“', Mz- = 3.5 GeV. Again 
usin g eq uation (fTOl i. The horizontal lines denote the range of Q/z^ measured by Planck Collabora¬ 
tion 113611 . 
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Figure 6; Relationship between the parameters g,j and mjj giving the correct DM relic abundance 
within 3(T. 


The y-ray flux summed over all possible annihilation ehannels leading to y-ray 
produetion is given by 


d^DM 

dE^ 


X 


1 

Anm^ 


'Eiio-iV) 


dN‘^ 

dEy 


1 



p^{r{s)\ds. 


( 20 ) 


where denotes the mass of the DM eandidate and, (crv) denotes the thermal 
averaged annihilation eross seetion. Sinee in our model m,, > Mz’, the dominant 
annihilation ehannel, as already mentioned, will be 77 + ^ > Z' + Z' This is then 
followed by Z' deeay. Thus, in our ease, the resulting y-ray speetrum is then the 
sum of all Z' deeay ehannels weighted by their branehing ratios and eonverted 
from hadrons and muons to y-rays. The integrals in equation (| 7 ]) are integrals over 
the line of sight to the target and averaged over the solid angle of the region of 
interest, AD. 

Various profiles of the radial distribution of dark matter in the galaxy have 
been put forth, see e.g., [37], and lli^ . We have used mierOMEGAs paekage 
and adopted the profile given and diseussed Refs, in and [39] to perform our 
ealeulations. We thus take a DM halo density profile p(r) normalized to a value of 
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Po at the solar galactic radius of 0.3 GeV/cm^. The function p(r)] is given by 1138h 


Pi.r) = Po (— 


Rc.\Ur' 


+ K 


ra + fO 


P-y 


( 21 ) 


We take Rq, the solar galactic radius, to be 8.5 kpc. The other parameters in 
equation (|2TI) are taken tobeQf = l,jS = 3,y = 1. 

We explore the production of WIMP annihilation y-rays obtained with our se¬ 
cluded DM model by following the analysis given in Refs. [[l^.[ll]l. In Ref. [[ill] it 
is concluded that the observed y-ray spectrum is best fit by WIMPS with a mass ~ 
20 - 50 GeV that annihilate to quarks with a cross section < (TamV > = 
cm^s“'. Motivated by that work, we consider here two specific cases for our 
seculded WIMP models, viz., WIMPS with masses of 20 GeV and 32 GeV. With 
the parameter choices for these models as given and discussed in Section [6l our 
results are compatible with the observed DM density as parametrized by the value 
Q.h^ ^ 0.119, and also with the value for {(Tannv) required to explain the proposed 
y-ray excess from DM annihilation. We now proceed to demonstrate this more 
explicitly for the two WIMP mass models that we consider. 


Model 1 : In this first example, we take m,, = 20 GeV, with the parameters as shown 
in the caption of Figures 2 and 3 chosen to fit the DM svalue for « 0.119 and 
the value for {cTannV) - 1.8 x 10“^® cm^s. The lifetime of the light mediator, Z', 
must be less than one second in order to guarantee that the Z' decays before the 
epoch of big bang nucleosynthesis ll22ll . We have chosen My = 0.5 GeV, which 
implies that Mz' - 353 MeV. Given this Mz' mass and all the couplings of Z' 
with fermions (see appendix), we then find the Z' decay width, Yz' - 1.76 x 10“^^ 
GeV. The branching ratios for the Z' decay channels are found to be 33.9% into 
uu, 25.4% into 24% in 8.48% into dd and 8.18% into ss. 

In this case, the annihilation y-rays primarily result from the decay of the Z' 
into light quarks uu,dd, ss, which accounts for 50.6% of the decay width. It is 
shown in Ref. lllin that the annihilations of DM particles of mass between ~ 18 
and ~ 26 GeV into light quarks may significantly account for the ~ 2 - 3 GeV 
excess. 

Choosing the set of parameters given by g^j = 0.088, m,, = 20 GeV, M^/m,, —> 
0, gvB = 8 X 10“^, and taking into account the constraint for tan 26a, designated by 
t 2 a in equation (ITOl) . we obtain Sa - -3.8 x 10^^. In this case, the spin independent 
cross section of rj with the nuclei is « 2.31 x 10“^^ cm^ and ~ 5.42 x 10“^^ 
cm^. The value of is of the same order of magnitude as that of the upper lim¬ 
its obtained by the both the LUX (Large Underground Xenon) experiment and the 
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XENON 100 experiment for the WIMP mass range that we eonsider here iUn 
our ealeulations, we obtained our parameters using mieromegas paekage mM. 
We note that and are dependent on gvB parameter, whieh in praetiee eon- 
neets the dark and visible seetors of the model (c/, Model 2 below) 


Model 2: In this model, we take = 32 GeV. We ehoose the interaetion 
strength = 1.12 x 10"', we take gvB = 1 x 10"^, and taking into aeeount the 
eonstraint for t 2 a in equation (fTOl) . we obtain Sa - -4.7 x 10“^. Then we inerease 
the mass of the Z' mediator to 3.5 GeV. The dominant interaction channel is still 
7] + 1] ^ Z' + Z'. In this case, we find (o-annv) - 1.88 x 10“^^ cm^s"' and we fit the 
DM ^ 0.119 from the results shown in Figures |4] and [51 Our values for the 
spin independent cross sections of rj with the nuclei are ^ 6.04 x 10"^® and 
o-'Ij - 1.42x 10"^^ cm^. For Mz' = 3.5 GeV, the additional cc decay channel opens 
up, being now kinematically allowed. The Z' decay width, Tz', will therefore be 
larger than that for Model 1, in this case Tz' - 3.74 x 10"'"' GeV. The branching 
ratios for the Z' decay are found to be: 25% in uu, 24.9% in cc, 18.8% into 
e'^c", 18.8% into 6.28% into dd, 6.28% into si. Thus, in addition to the 
other channels of Model 1, the cc now also significantly contributes to 
yy production. For a large, cc decay channel it is shown in Ref. [fnl] that the 
annihilations of DM particles of mass between ~ 28 GeV and ~ 36 GeV into cc 
channels may significantly account for the ~ 2 - 3 GeV excess. 

The expression for Z' decay width into fermion-antifermion pairs, //, is given 
by 


r 

^ z'^// 


= z 


1_L 

c\ lln M\, 




4m? X 


- Ami) + + 2 m?)], 


( 22 ) 


where the sum runs over quarks and leptons species kinematically allowed. Nc 
is the color number, m, denotes the mass of the fermion and y denotes ax¬ 
ial/vectorial couplings of Z' to fermions. Therefore, in the model considered here, 
the Z' will decay more into up-type quarks than down-type quarks. Even if we 
increase Mz’ so that Mz' > 2m/,, the bb channel will still not dominate. We could 
open this channel kinematically, but this would require a larger mass for rj in order 
to keep the primary annihilation channel as to be 77 - 1-77 ^ Z' -l-Z' in the mass range 

» Mz’. 

If, for example, we take the same set of parameters for Model 2, only increas- 
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Figure 7: Predictions of the excess y-ray flux from WIMP annihilation for our two models as 
described in the text (green: Model 1, orange: Model 2). The upper limits of the energy spectrum 
are determined by the WIMP masses since annihilation occurs near rest. 

ing Mz’ ~ 9.28 GeV, so that Z' —> bb is now opened, than the branehing ratios 
of Z' would be 20.1% in uu, 20.1% in cc, 15% into e'^e~, 15% into 14.9% 
in 5.15% into dd, 5.15% into ss and only 4.57% into bb. One of the main 
phenemologieal differenees between this model and the other ones presented in 
the literature is that, as ean be seen in this example, Z' ^ bb will not the most 
probable deeay ehannel. The detailed expression for y ean be found in the 
appendix. 

Considering the values for ehosen for model 2, for instanee, if Mz' > 2mr^, 
the Z' deeays 100% into 7777 . However, smaller values for ean suppress this 
deeay and ehange this ratio. The expression for Z' —> rjfj is given by 



Figure 0 shows the y-ray flux as a funetion of the energy of the photon for 
the speeifie models that we eonsidered. As with most previous work, we negleet 
a possible seeondary eontribution from Compton seattering of the eleetrons and 
positrons produeed in the annihilations, although it has been suggested that this 
proeess may eontribute to the y-ray exeess ll^ . 

8. Conclusions 

We have proposed a model in whieh a Proea-Wentzel (PW) field eonsisting 
of ~ MeV-multi-GeV mass partieles are mediators of seeluded fermionie DM 
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interactions. From the low-energy phenomenology it is known that electroweak 
precision data such as the mass of W boson, the decay width of Z boson, and 
some asymmetries, can constrain such a model. For our secluded DM model, 

the constraints on -^|1 - /M||, where ^ = gvB/ - gvB’ satis¬ 

fied. Indeed, they are much smaller than 0(10^^) 112311 . Also, the mediator of DM 


interactions will not induce flavor-changing neutral current processes. 

One important difference between our model and other ones presented in the 
literature is the absence of an extra singlet scalar. Thus, there is no equivalent to 
a ’’Higgs scalar” in the PW Lagrangian. This distinguishing characteristic may 
become important if no /z ^ invisible width is observed at LHC; such an result 
could be a clear signature of the model. 

The model that we have explored has characteristics in common with the 


model presented in 01911 . Our results agree with theirs in the limit Mz' > 0. As 


those authors point out, there are interesting and testable signatures for secluded 
DM models. We now summarize various empirical implications of our secluded 
DM model, based on the results of our calculations. 

8.1. Indirect tests: Gamma-ray Excess 

We have here explored one of the testable signatures of our secluded WIMP 
models in quantitative detail. We have considered the y-rays that would be pro¬ 
duced as a result of cosmic annihilations of secluded WIMPs of mass of ~ 20 and 
~ 32 GeV. In particular, we have shown that the secluded DM model proposed 
here can potentially explain both the flux and spectrum of an apparent 2-3 GeV 
energy y-ray excess in the galactic center region, this being an excess over that 
expected by taking account of other galactic y-ray production processes. Such an 
excess has been inferred from an analysis of Fermi-hAI Gamma-ray Space Tele¬ 
scope data Our results are also consistent with a putative weak ~(2-3) cr 
y-ray signal claimed for the dwarf galaxy Reticulum 2 [|43ll and with the conser¬ 
vative constraints derived from Fermi Gamma-Ray Space Telescope observations 
in Refs. [441. 

Unfortunately, owing to the uncertainties in the empirical derivation of a y- 
ray excess, it is difficult to distinguish between the y-ray spectra obtained for our 
WIMP models and those that generally follow the analysis of Ref. [Q . However, 
we note that most other DM interpretations of the galactic center excess stress the 
bb channel, whereas in our Model 1 light quarks and leptons dominate in produc¬ 
ing the DM annihilation y-rays while in our Model 2 the cc channels dominate 
over the bb channel. This is a predicted phenomenological difference between 
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our model and other models. At present, owing to the observational uneertain- 
ties in determining the DM ’’signal” over the other proeesses eontributing to the 
y-ray ’’baekground” in the direetion of the galaetie eenter, a definitive test of this 
differenee is diffieult. 

We note that the DM interpretation of the ’’exeess y-rays from the galaetie 
eenter region is subjeet to some eaveats. We again note that there are other inter¬ 
pretations of either all or part of this exeess [Il3n . Also, using a new ealibration 
ealled Pass-8, the Fermi eollaboration has very reeently argued for upper limits on 
the eross seetion for the produetion of y-rays from DM in dwarf galaxies that is 
in tension with the DM interpretation of the galaetie eenter exeess [l45n . However, 
this involves modeling of the distribution of the square of the density of annihi¬ 
lating partieles along the line-of-sight in these galaxies [|4^ . now referred to as 
the ”J-faetor” [|47ll . These models, in turn, rely on veloeity dispersion measure¬ 
ments [|48ll that do not distinguish between the gravitational elfeets of DM and 
faint stars. 


8.2. Indirect tests: Cosmic Ray Antiprotons 

The AMS-02 experiment has obtained detailed data on the eosmie ray antipro¬ 
ton speetrum in the vieinity of the Earth ll^ . The implieations of these results 
for DM annihilation models involve various uneertainties sueh as those of galae¬ 
tie eosmie-ray propagation, eross seetions for interstellar p produetion by eosmie 
rays, and solar modulation. The authors of Ref. ll^ have eoneluded that these 
data provide no unambiguous evidenee for a signifieant exeess above that ex- 
peeted from eosmie ray interaetions with interstellar gas. Should p produetion by 
DM annihilation be absent, we wish to point out that this would favor our Model 1 
with a value of Mz> below the kinematie threshold for p produetion from Z' deeay. 
On the other hand, p produetion is kinematieally allowed in our Model 2. Thus, 
should future observations definitively rule out an observable p eomponent from 
DM annihilation, this would support a seeluded WIMP model for explaining the 
galaetie eenter y-ray exeess. 


8.3. Direct Tests 

For both of the seeluded WIMP model parameters that we eonsider here, our 
ealeulated values for the eross seetions for WIMP-proton seattering are eonsistent 
with the present experimental upper limits obtained by the XENON 100 and EUX 
experiments The WIMP-neutron eross seetions are suppressed. It is antiei- 
pated that by 2020 liquid xenon deteetors will have the eapability to measure spin 
independent eross seetions as low as OCIO""^^) em^ [15 ill . Should future laboratory 
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results yield a very small eonstraint on the WIMP elastie seattering eross seetion, 
the DM annihilation hypothesis for explaining the y-ray exeess from the galactie 
eenter region would then favor a seeluded WIMP model for the dark matter. 
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Appendix A. Interactions and vertices of the model 

Couplings of Z and Z' to SM fermions: 


l//lf/Z . Cq.(1 

Cw 

_ ig 

if/l^Z . Caita 

Cw 


Tl- 




n- 


(1 

>H/( 






(ta + Sw^) ' 

Triple gauge bosons eouplings: Comparing to the SM eouplings denoted by 
R, they will be: 


Higgs Couplings: 


Raw+w- = C 
Rzw+w- = Ca, 
Rz’W+W- = -Sa- 


nif 

hff : 

Vh 

,2 


m 


hWW : 2i^, 
Vh 


mx 


hZZ : 2i—(Ca - ^SwSaf', 
Vh 


mx 


hZ'Z' * 2/- {Sa "I" ^S\i/C(y) , 

Vh 


mx 


hZZ' : 2i —[(c„ - -t Sw^cJ]. 

Vh 


(A.2) 

(A.3) 
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Coupling of Dirac Fermion and Z and Z' gauge bosons: 


?7/?7Z: igr^ =Sa, (A.4) 

rifr]Z' : igr, ^ =c„. (A.5) 

The interactions of Z' in Eq. (lA.ll) are written in a simplified form. In faet, in 
order to calculate the Z' decay width we have used eouplings in the form of 


-Cnc = 'y.Wi'fig'v - - fAy^)4'iZ'f,], (A.6) 


where gy, g\ denote respeetively the veetorial and veetor-axial coupling of 
Z boson with fermions and fy, denote these couplings but now for Z' gauge 
boson. The expressions whieh relate vectorial and vector axial couplings to left 


and right-handed eouplings are given in Eq. Appendix A 


gA 

gv 

fk 


1 

2^gL ~ gR)’ 

(A.7) 

1 / 

2(gL + gR)’ 

(A.8) 

2^fL ~ /«)’ 

(A.9) 

5(/i + /«)- 

(A.IO) 
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The detailed left and right-handed couplings of Z' to fermions are given below. 


/“ = cMa+^S„) ^ (A. 11) 

\ta + ^Sw) 

Ir = Ca(ta+^Sw) ’ (A.12) 

yta+^Sw) 

ft = (A.13) 

\^a ' S 

Ir ~ C(y{ta ft I i: \ ^wQd 5 (A.14) 

{ta + ^Sw) 

fl = cMa+fsw)\Tl- ^lt^^'/'^\ lQy (A.15) 

Va + ^Sw) 

Ir = Ca(ta+fSw) ---, (A.16) 

{ta + ^Sw) 

fl = cMa+fsw)Tl (A. 17) 

fl = 0, (A.18) 

where Qi denotes the charge of fermion, = -1/2, TJ = 1/2. 

For completeness we write here the right and left-handed couplings of Z to 
fermions. 

8l = [C.(l - Swta^)] \t; - ~ . (A.19) 

gl = [Caa-Swta^)]\-jf^^^^slQl (A.20) 

gi = [c.(l - 5 wt.^)] \tI - ~ , (A.21) 

(i S\YtQ;g) 

gl = [Ca(l-swta^)]\-jf^^^^slQd], (A.22) 

gl = [C.(l - Swta^)] It; - ~ , (A.23) 

(1 ^W^ab) 

gl = [Ca(l - Swta^)] - n , (A.24) 

gl = (c^-Swtam;, (A.25) 

gl = 0. (A.26) 
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